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We consider the optimal conduction path of the one-dimensional variable-range hopping problem. 
We describe a hierarchical procedure for constructing the path which is in excellent agreement 
with numerical results obtained from a percolation approach. The advantage of the hierarchical 
construction is that it is easier to analyse. We show that the distribution of hopping lengths is 
well approximated by a model for the repeated breaking of a stick at its weakest point, until the 
fragments are too strong to be broken. 



PACS numbers: 05.40.-a,05.60.-k, 72.20.Ee,71.55.Jv 

1. Introduction. The concept of variable-range hop- 
ping was introduced by Mott [l| to explain the empir- 
ically observed temperature dependence of the electri- 
cal conductivity in disordered semiconductors at very 
low temperatures. Mott argued that the conduction is 
determined by a competition between transitions with 
large matrix elements and transitions requiring a small 
activation energy: large matrix elements favour short- 
range hopping, but there are unlikely to be energetically 
favoured transitions at short ranges. Maximisation of 
the transition rate suggests that the conductivity should 
have a temperature dependence of the form [l[ 



cr(T) - Ae 



-(T /T) 1/(d+1) 



(1) 



It is very difficult to produce accurate quantitative results 
from Mott's heuristic arguments. Ambegoakar, Halperin 
and Langer Q pointed out that the resistance of a sam- 
ple is determined predominantly by the highest resistance 
of the 'optimal conduction path', and demonstrated how 
Mott's picture is related to percolation. They showed 
that the conductance is of the form of (Q]) for d > 1, 
and could relate T to the geometry of percolating con- 
ducting clusters, which can be determined numerically. 
Later Kurkijarvi Q pointed out that in one dimension 
the mean value of the resistance is determined by acti- 
vation over the most unfavourable paths, and that for 
(one-dimensional) wires, the mean resistance is of Arrhe- 
nius form at very low temperatures. His argument was 
improved upon by Raikh and Ruzin Q who calculated 
the exponential term in the mean resistance accurately. 
Thus the only precise results concern the weakest link in 
the conduction path, despite the fact that this problem 
has been investigated for four decades. It is desirable 
to gain a more complete picture of the conduction path. 
In this paper we consider a quantitative model for the 
distribution of hopping lengths in one dimension. 

Efros and Shklovskii [H have suggested that electronic 
correlation effects could modify the predictions of Mott's 
variable-range hopping picture. Mott's approach is usu- 



ally discussed in the context of a degenerate electron gas, 
in which case correlation effects may be significant. How- 
ever we consider a slightly modified problem applicable 
to non-degenerate systems, thus avoiding possible com- 
plications due to correlation effects. 

2. Results. In the limit T — ► the one-dimensional, 
non-degenerate variable-range hopping model is mapped 
to a resistor-network model (Fig. [I} with conductances 
T mn . As T — ► 0, conduction is confined to a single most 
favourable path. Our aim is to determine the properties 
of this 'optimal conduction path'. We characterise this 
path by the distribution of its hop lengths. 

We introduce a hierarchical approximation scheme for 
the optimal conduction path. An example is shown in 
Fig. In this example, the hierarchical construction 
yields exactly the same path as a numerical percolation 
approach (described below). For the parameters in Fig. 
[2k, the hierarchical construction yields the correct hops 
in 98.2% of all cases. The corresponding distribution of 
hop lengths is shown in Fig. 03, in excellent agreement 
with results from a percolation approach. The advantage 
of the hierarchical procedure is that it is susceptible to 
geometrical and statistical analysis. We show that the 
problem can be analysed in terms of repeatedly breaking 
a stick at its weakest point until the fragments cannot 
be broken down further. The distribution of fragment 
lengths is found to be in very good agreement with the 
distribution of hopping lengths (Fig. [2J:) . 

By microscopic probes (such as those used in tunnel- 
ing microscopy) it is becoming possible to determine the 
paths involved in carrying the current in microscopic sys- 
tems. The distribution of hopping lengths may thus be- 
come accessible to experiments. 

The remainder of this letter is organised as fol- 
lows. First we consider a resistor-network model for the 
one-dimensional, non-degenerate variable-range hopping 
problem. We then describe our percolation approach to 
determining the optimal conduction path (an extension 
of the approach employed in 0). Our results serve as 
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FIG. 1: One-dimensional resistor-network model for the one- 
dimensional variable- range hopping problem described by 
equation The conductances F mn are given by ((3j), E n 

are random on-site energies. 

a benchmark for the hierarchical approach which is de- 
scribed next. The remaining sections show how to com- 
pute the distribution of hop lengths in terms of a stick- 
breaking process. 

3. Resistor-network model. We consider an array 
of sites with energies E n , independent identically dis- 
tributed random variables in [0,1]. We assume that the 
occupation probabilities P n satisfy a rate equation (ne- 
glecting quantum-mechanical interference effects) 

Pfl ^ ^ RmnPm Pn ^ ^ Pnni • (2) 
m m 

The transition rate from state n to state m is 
called R nm ] the principle of detailed-balance implies 
Rnm = R mn cxp[-(E m - E n )/(k B T)}. We specify the 
rates by assuming that downward transitions {E m < E n ) 
occur at a rate R mn = el"~ m l for e > 0. 

The master equation ([2]) is inconvenient in its orig- 
inal form, because its equilibrium distribution P„ = 
Z" 1 exp[— E n /(k-B,T)] is non-uniform. We therefore 
transform it to an equation evolving towards a uniform 
equilibrium density, and then describe deviations from 
this equilibrium in terms of equations analogous to Kir- 
choff's laws for a resistor network, with effective con- 
ductances T mn between lattice sites m and n. We then 
expect (following the line of argument employed in 
that the diffusion rate is determined by a resistor network 
(Fig. [l| with 'conductances' T rnn = Pm^ Rnm'- 

T„ m = e \n-m\ e -max(E m ,E n )/(k B T) = e ~A mn /t _ (3) 

The last equality defines 'distances' A m „ = at\n — m\ + 
max(i? m , E n ), with a — — loge, and t = k^T. 

4. Percolation approach. In the limit as t — > 0, ratios 
between the conductances T mn become very large. In 
this limit the fraction of the current in any bond is almost 
certainly very close to zero ('inactive bond'), or else very 
close to unity ('active bond'). In [3] it was pointed out 
that the resistance of the network is almost surely well 
approximated by the resistance of the active bond with 
the highest resistance. This bond can be identified by 
mapping to a percolation problem: for a given choice 
of threshold conductance To, bonds with conductances 
r m „ <To are removed from the network. Decreasing To, 
we test for connectivity of the network. As To decreases 
below the conductance of a 'critical' bond, (n, m)* (with 




FIG. 2: a One realisation of the hierarchical process described 
(a = 0.1, t = 0.001, Af = 2000). When segments break they 
do so mostly by one break (thick black vertical lines). In one 
instance a segment breaks by two simultaneous breaks (red 
lines) into three pieces. The process terminates in generation 
8 when all segments are unbreakable. For the particular reali- 
sation shown, the percolation approach gives exactly the same 
segments, b Distribution of hop lengths of the optimal con- 
duction path (■) determined from a percolation approach. 
Also shown are approximations employing the hierarchical 
process, allowing for one (o) and two (0) breaks per segment. 
Allowing for more than two simultaneous breaks per interval 
does not change the distribution significantly. The parameters 
were a = 0.1, t = 0.002, J\f = 5000. c Stick-breaking approxi- 
mation for Phop (AO allowing for one break per segment (with 
the same parameters as in b). We plot equation (111[) . black 
line, and on two iterations of l)10p. red line. Further itera- 
tions do not change the red line noticeably. For comparison 
we show data from the hierarchical scheme allowing one break 
per segment (o), as in b. 



conductance r* TO ), we find that the network becomes 
connected. The conductance of the network is expected 
to be well approximated by T* mn . Now we repeat this 
procedure for the sub-networks to the left and to the 
right of the critical bond, determining the critical bonds 
of the sub- networks. Repeating this procedure results in 
identifying critical bonds in progressively shorter sections 
of the network. The process terminates when the critical 
bond for a sub-network is the one bond that connects its 
ends directly. In this way the optimal conduction path 
is found. This approach is accurate when t — > 0, but the 
search over all possible paths is very time consuming, and 
not susceptible to statistical analysis. Our objective was 
to find a way of determining the optimal conduction path 
by repeated application of simple rules - which can be 
analysed statistically. 

5. Hierarchical model. Our approach consists of 
breaking down the network by proposing a succession 
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of current-carrying paths with smaller and smaller resis- 
tances. At each stage we consider a segment (to, n) of 
the network, between nodes to and n. We ask whether 
the conductance of the direct connection T mn is less than 
that of an indirect path, which goes through an interme- 
diate node j. If the resistance T^j+Tj* is less than r~^, 
then the network is broken at node j, and we consider 
the sub-networks (m, j) and (j,n). If there is more than 
one break which lowers the resistance, we choose the one 
which gives the lowest resistance. If there is no single 
'break' lowering the resistance, we try inserting two si- 
multaneous breaks, then three, etc. If there is more than 
one possible way to insert the breaks, the choice giving 
the lowest resistance is used. If the resistance cannot be 
reduced below that of the direct connection by consider- 
ing any sequence of nodes connected in series, then the 
sub- network is determined to be 'unbreakable'. 

There is no guarantee that the process will produce 
the optimal conduction path, only that the resistance will 
decrease with each refinement. As Fig. [2]shows, however, 
this procedure provides an excellent approximation to the 
optimal conduction path. 

6. Stick-breaking model. The hierarchical method has 
the advantage that it is easier to analyse than the perco- 
lation approach. In the remainder we illustrate this by 
showing that the simplest hierarchical model (allowing 
only one break per segment at a time) can be approxi- 
mated by a 'stick-breaking model'. 

It is found that when an interval can be broken at one 
point to reduce the resistance, the break point usually 
corresponds to the site within the interval with the small- 
est energy Ej. This observation leads to the following 
simplified model for the hierarchical process. When con- 
sidering an interval between sites n and to (with n > to) , 
we determine whether the resistance is reduced by in- 
serting a single break. If so, we determine the smallest 
energy Ej with to < j < n, and break the interval there. 
The process of sub-division is continued until all of the 
intervals are determined to be unbreakable. 

In order to model this process we introduce two prob- 
abilities. Let Pbr(AQ be the probability that a fragment 
of length N elements is breakable. If it is breakable, it 
will break at a position M elements from the end with 
probability p(M,N). Starting from and infinitely long 
stick and iterating this process, we end up with a set 
of unbreakable fragments with a probability Ph& g (N) for 
having length N. In the following, we discuss how Pt, T (N) 
and p(M, N) respectively are chosen to make this stick- 
breaking model correspond to the simplified hierarchical 
model of variable range hopping. We then discuss how 
the stick-breaking model is solved. Figure [2h shows that 
the resulting distribution of fragment lengths, Pi ia ,g(N), 
is a good approximation to the distribution of hopping 
lengths, P hop (N). 

7. Probability for a segment to break. The set of 
energies can be represented by a set of points in the j, Ej 
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FIG. 3: Shows points in (j, Ej) plane (j = 0, 1, 2, . . . and Ej 
uniformly distributed on [0, 1]). The probability that the seg- 
ment (m, n) is unbreakable (with a single break) is equal to the 
probability that the 'roof of the house' (hashed) of area A' is 
empty. The dashed horizontal line represents max{_E m , E n }, 
the second dashed line Ej = max{_B m , E n } + at (j — m). 

plane, which are randomly distributed in energy with the 
density of points in the plane being unity. 

In the limit of t — > 0, the resistance of an indirect 
path through one intermediate node is approximated by 
the larger resistance of the two links, so that we write 
r ny + r 7« ~ exp[max(A mj -, Aj n )/t]. So the probability 
for a segment [to, n] to be unbreakable is the probability 
that A n j > A nm and Aj m > A nm for any choice of j 
(with to < j < n). These two conditions imply that the 
regions below two lines in the (j, Ej ) plane are empty, in 
the interval [to, n]. Thus the probability that the interval 
[to, n] is unbreakable (with a single break) is equal to the 
probability of finding the 'house' in Fig. [3J empty. 

The probability of finding a region in the (j, Ej ) plane 
of area A to be empty of points is exp(— A). At first 
sight it seems as if we should take A to be the area of the 
'house' in figure [3J and the probability that a segment is 
breakable would be Pbr = 1 — e~ A . We do however have 
some prior knowledge about the energies: because we 
always break at the minimum energy, the none of the en- 
ergies Ej in the interval [to, n] is less than max(-E„, E m ). 
Thus it is the area of the 'roof of the house (hashed in 
Fig. [3J which is relevant, and the probability for this 
segment being breakable is therefore 

P hl (n - m) = 1 - exp(-A') , A' = at(n - to) 2 /4 . (4) 

Thus the probability for a segment being breakable in our 
construction is, in fact, solely dependent on its length, 
and given by (@|. 

8. Distribution of break points. Next we consider the 
probability p(M, N) that an interval of length N will 
break at a position M units from one of its ends. Be- 
cause the energies Ej are independently and identically 
distributed, we might expect that the position of mini- 
mal Ej is uniformaly distributed, that is p(M, N) should 
be independent of M. This is incorrect, however, be- 
cause we are considering the distribution of energies in 
an interval which is selected to be breakable. 

The probability p(M, N) for the position M of mini- 
mal Ej in an interval of length N is symmetric around 
the midpoint of the segment. Let x S [0, 1] be a vari- 
able measuring distance of the break from the end of 
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the interval, setting x = 2M/(N - 1) if M < N/2 and 
x = 2(N - M)/(N - 1) if N/2 < M < N. We consider 
the continuous probability density p(x) for x, which is 
related to p(M, N) by: 



p(M,N) 



1 J p(^) ifO<M<AT/2 



AT- 1 



p(^f) N/2<M<N. 



(5) 



The density p(x) is calculated by geometrical analysis 
of Fig. [3] We rescale the energy Ej to a new variable 
y = (N — l)[Ej — max(E n , E m )] so that the density of 
points in the (x, y) plane is unity. In these new coor- 
dinates the line determining the boundary of the region 
which must be occupied to ensure breakability is y = vx 
with v = atN 2 /2. Instead of calculating p(x) we com- 
pute the complementary probability density q(x), that is 
the probability density for the position of the minimum 
energy in an unbreakable segment. The probability den- 
sity for the location of minima irrespective of breakability 
is uniform (and by normalisation equal to unity), so that 
1 = Pbrp(x) + (1 - Pbr)q(x), or: 

p(x) = [1 - exp(-v/2)q(x)]/[l - exp(-i//2)] . (6) 

It is easier to compute q(x) because the conditional 
knowledge about unbreakable (at one break) intervals is 
simpler. These intervals have no points in the (x, y) plane 
inside the 'forbidden' triangle < y<vx, 0<x<l. Thus 
q(x) is the distribution of x for the point with minimal 
y from a random, unit density scatter in the 'allowed' 
region y > vx, < x < 1. Let F(y ) be the probability 
that the point with lowest value of y is above y$. This 
satisfies F(0) = 1 and dF/dy = -dA/dyF where A(y) 
is the area of the allowed region below y. We have 



q(x) 



dF 
dy 



X[o,x(v)]( x ) 



x(y) 



(7) 



where X[a,b] ( x ) is the characteristic function (equal to 
unity on the interval [a, b] and zero elsewhere) and x(y) — 
dA/dy is the width of the allowed region at y. Us- 
ing the differential equation for F(y) we obtain q(x) = 
J™ dy F(y). We find F(y) = e -v 2 / 2v for < y < v and 
F(y) — e~ y+u l 2 for y>v and finally obtain: 



q(x) 



-u/2_ 



erf 



, erf ( x , , 
2 V V 2 



(8) 



Eqs. ([5][5]) show that segments are more likely to break 
in the centre, as expected. This bias is the stronger the 
shorter the segments are. For long segments (v — ► oo) 
we obtain p(M,N) w 1/(N - 1). We find that a good 
approximation to q{x) is 



q(x) w 1 +v(l/2 - x) . 



(9) 



9. Self-consistent solution of the stick-breaking pro- 
cess. We now determine the probability distribution 



Pfrag(N) of the lengths of the unbreakable fragments, by 
extending a calculation for the distribution of strength of 
a repeatedly broken random chain, discussed in [2]. We 
start from a very long stick, and after T steps of this 
process, we have Wd(N,T) segments of length N which 
have been determined to be unbreakable, and W\j(N, T) 
segments of length N which have not yet been tested. As- 
suming that these numbers are so large that it is sufficient 
to consider expectation values, we have the recursion 
W V (N,T+1) = 2j2^ I=N+1 p(N,M)PUM)W v (M,T) 
andW D (N,T+l) = [l-P bt (N)]W v (N,T). Rather than 
following this iteration for a single stick being broken, we 
consider a steady state of W\j (N, T) , with destruction of 
one additional stick being initiated at each step: 



W(N) = 2 P( N > M)P hl (M)W(M) . 



(10) 



M=N+1 



The corresponding probability distribution of fragment 
lengths (or hop lengths) is Pha.g(N) = Z~ x \l — 
Pbi(N)]W(N) where Z is a normalisation factor. 

When p(N,M) = 1/(M - 1), equation HUJ) is solv- 
able. Replacing sums by integrals, we obtain Wq(N) oc 
iV -2 e ~ El ( atN / 4 ) j where Ei is the exponential integral. 
When p(M,N) is given by eqs. ©, ©, and ©, we 
solve (|10[) by iteration starting with Wo. Usually a few 
iterations give an accurate solution. The first iteration 
gives [using eq. © to approximate ©] 



ffr ag (iV) cx e~ v ' 2 [h{ V )+I 2 { V )] 



h(v) 



d/z 
772 

Av A 4 



^(e^-l + ^-2V^)e 



2 



(11) 

-Ei( M /2)-/i/2 
-EiQV2)-M/2 



The distribution P[ rag (N) given by (jlip is shown in Fig. 

Results obtained from higher iterations of (|10[) using 
([5]) are also shown. They converge rapidly and provide a 
strikingly accurate approximation to Phop(^0- 
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